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ELECTROSTATIC  WAVE  OBSERVATION  DURING  A 
SPACE  SIMULATION  BEAM-PLASMA-DISCHARGE 


I.  INTRODUCTION 

The  Beam-Plasma-Discharge  (BPD)  instability  (Getty  and  Smullin,,  1963)  has 
re-emerged  in  recent  years  (Smullin,  1981;  Bernstein  et  al.,  1978)  as  a  funda¬ 
mental  process  to  be  considered  by  designers  of  space  vehicle-borne  electron 
accelerator  experiments.  Because  of  concerns  over  vehicle  neutralization 
during  these  experiments  and  an  interest  in  further  understanding  of  the  BPD 
process  (Cambou,  1975,  1978)  a  series  of  space  simulation  investigations  of 
the  BPD  phenomenon  was  undertaken  in  the  large  vacuum  chamber  facility  at 
Johnson  Space  Center  (Bernstein  1978,  1979;  Szuszczewicz  1979,  1982a; 
Papadopoulos,  1982;  Anderson,  1981).  Bernstein,  et  al.  1979  were  able  to 
conclude  that  the  phenomenon  is  dependent  upon  a  critical  electron  beam 
current,  Ig  ,  for  ignition  and  that 

1°  *  Vb3/2  B”*7  p“1l_1  (1) 

where  V0  i3  the  beam  energy  in  volts,  B  the  magnetic  field  in  gauss,  P  the 
ambient  neutral  particle  pressure  in  torr,  and  L  the  beam  length  in  meters. 

The  BPD  is  characterized  by  a  transition  from  single  particle  to  collec¬ 
tive  processes.  Further  characteristics  of  BPD  include  enhanced  electron 
densities  and  spatial  broadening  of  the  beam.  In  addition,  there  are  a  number 
of  plasma  mode  oscillations  excited  which  extend  from  plasma  waves  at  frequen¬ 
cies  near  the  local  plasma  frequency  down  to  electron  cyclotron  and  extremely 
low  frequency  (ELF)  electrostatic  oscillations.  In  a  previous  paper 
(Szuszczewicz  et  al.,  1982a)  we  established  a  density  dependent  criterion  for 
BPD  ignition.  In  the  present  paper  we  examine  ELF  density  oscillations 
along  with  plasma  potential  and  frequency  measurements  and  discuss  a  model  of 
azimuthal  wave  propagation  which  is  consistent  with  our  observations. 
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Low  frequency  drift  wave  instabilities  have  long  been  associated  with 
energy  transport  and  anomalous  diffusion  in  turbulent  plasma  regimes  and 
have  hence  been  an  important  topic  in  studies  of  plasma  confinement 
(Kadomtsev,  1965).  In  low-8  space  applications  the  sharp  density  gradients 
encountered  in  the  equatorial  ionospheric  studies  of  spread-F  (Ossakow,  1981; 
Huba  and  Ossakow,  1979;  Singh  and  Szuszezewicz,  1984;  Kelley,  1982;  Gary, 

1983;  Szuszezewicz,  et  al.,  1982b)  are  often  sufficient  to  produce  low 
frequency  drift  instabilities.  In  the  laboratory  environment  described  here, 
earlier  investigations  of  the  BPD  phenomenon  (Szuszezewicz,  et  al.,  1979) 
yielded  estimates  of  perpendicular  diffusion  coefficients  orders  of  magnitude 
higher  than  predicted  for  classical  collisional  diffusion.  Those  results  were 
consistent  with  a  Bohm-like  diffusion  process  created  by  wave  turbulence.  The 
current  experimental  configuration,  described  below,  is  an  outgrowth  of  these 
results  and  a  recognition  of  the  desirability  of  further  quantifying  the  wave 


measurements. 


II.  EXPERIMENTAL  CONFIGURATION 


The  experiment  was  conducted  in  the  large  vacuum  chamber  at  the  Johnson 
Space  Center  with  a  configuration  similar  to  earlier  investigations  of 
Bernstein  et  al.,  1979.  The  primary  diagnostic  was  the  pulsed-plasma-probe, 
P^,  system  illustrated  in  Figures  1A  and  IB.  Figure  IA  is  a  schematic  view  of 
the  cylindrical  chamber  showing  the  position  of  the  electron  gun  and  collector 
plate.  The  magnetic  field  maintained  during  the  BPD  investigation  was  held 
constant  at  1.5  gauss  and  was  a  result  of  the  ambient  field  and  the  field 
produced  by  the  Helmholtz  coils  pictured  in  the  Figure.  Figure  IB  is  a  top 
view  of  the  chamber  showing  the  traversal  axis  of  the  P^  system  with  respect 
to  the  gun-collector  position  within  the  chamber.  A  complete  probe  system 
traversal  consisted  of  a  cycle  from  A  to  B  nd  back. 
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The  beam  was  generated  by  a  tungsten-cathode  Pierce-type  diode  gun  and 
was  aligned  as  closely  as  possible  to  the  magnetic  field. 

The  plasma  was  generated  by  beam  injection  into  the  ambient  neutral  gas 
(i.e.,  air)  at  pressures  near  10  ^  torr.  For  the  case  studied  here  the  beam 
energy  was  1.6  kev  and  beam  current  was  40  ma.  A  relative  electron  density 
profile  for  a  complete  traversal  during  BPD  is  shown  as  the  upper  plot  in 
Figure  2.  The  abscissa  is  time  measured  from  the  beginning  of  a  traversal  and 
the  vertical  lines  drawn  in  the  figure  identify  a  point  approximately  0.6m 
radially  outward  from  the  beam  center  and  effectively  demarcate  the  beam  core 
boundary  during  BPD.  In  addition  shown  in  Figures  2A-2C  are  radial  profiles 
of  electron  density,  electron  temperature  and  plasma  potential.  An 
exponential  least  squares  best  fit  to  each  data  set  is  also  shown  in  these 
Figures.  The  abscissa  is  radial  distance  with  respect  to  the  beam  center. 

III.  EXPERIMENTAL  TECHNIQUE  AND  RESULTS 

The  experimental  configuration  of  a  two-probe  system  in  which  the  probes 
are  one  meter  apart  along  the  radius  provides  a  2-point  measurement  for  phase 
delay  calculations  using  cross-correlation  signal  analysis  techniques  (Bendat 
and  Piersol,  1971). 

The  basic  operating  mode  of  the  two-probe  P^  system  provides  a  continuous 
two-point  determination  of  relative  density  variations  through  measurement  of 
electron  and  ion  saturation  currents  in  addition  to  the  simultaneous 
determination  of  absolute  densities,  temperatures  and  plasma  potentials 
through  characteristic  Langmuir  probe  response  to  an  applied  sweep  voltage. 
This  is  accomplished  through  a  pulsing  procedure  which  alternately  pulses  the 
probe  between  an  applied  sweep  voltage  and  a  fixed-bias  level  (electron 
saturation  for  the  E-Probe  and  ion  saturation  for  the  I-Probe).  In  this 
experiment  the  two-point  cross-correlation  measurements  were  done  using  E-  and 


I-probe  fixed  bias  levels.  Each  of  these  current  fluctuations  is  proportional 
to  density  and,  in  the  absence  of  potential  variations  or  ion  mass  changes, 
would  be  in  phase  in  a  homogeneous  neutral  plasma. 

The  basis  of  the  cross-correlation  technique  is  the  integral, 
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where  Igg(t)  is  electron-saturation  current  collected  by  the  E-probe  (at  a 
fixed  positive  bias)  and  Igj(t)  is  ion-saturation  current  collected  by  the 
I-probe  (at  fixed  negative  bias)  which  is  positioned  1  meter  radially  outward 
from  the  E-probe.  The  currents  are  detrended  using  a  least  squares  fit  to  the 
data  and  fluctuations  from  this  fit,  s*gg(i)*  were  normalized  by  the  fitted 
value. 

Figure  3  is  a  sample  of  a  .25  second  interval  of  data  taken  near 
r£  =  2.5  meters  from  beam  center  where  r£  indicates  the  E-probe  position 
relative  to  beam  center.  Illustrated  in  Figure  3  is  the  time  delay  tq 
(-/<  10  msec)  of  the  wave  train  passing  through  the  two-probe  system.  tq  is 
shown  more  clearly  in  the  plot  of  C(t)  in  Figure  4,  where  the  integration 
interval  was  1.2  seconds  and  the  sample  frequency  was  500  Hz.  Since  a  general 
waveform  such  as  C(r)  can  be  composed  of  a  number  of  prime  discrete 
frequencies,  it  is  necessary  to  decompose  the  waveform  in  order  to  uniquely 
associate  phase  shift  with  frequency.  The  Fourier  transform  of  C(t),  referred 
to  as  the  cross-spectral  density  function  (F  (i»)),  accomplishes  this 
decomposition  and  is  generally  a  complex  number, 

Fc(l.)  .  |Fc(i»>|.  ■1("'')IBE.IBt  (3) 
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In  this  form  (ut) 


is  the  phase  shift  of  the  frequency  component  w,  and 


IBE,IBI 

t  is  the  time  delay  of  propagation  through  the  two-probe  system.  The  cross 
spectral  density  in  Figure  5  shows  a  peak  near  45  Hz  and  a  phase  shift  near 
2.6  radians  consistent  with  the  time  delay  (r  )  in  Figure  4. 

It  is  appropriate  to  point  out  here  that  for  certain  non-local  propagat¬ 
ing  drift  modes  (Chen,  1967)  in  cylindrical  coordinates,  radially 
displaced  sensors  would  detect  a  phase  difference  in  density  oscillations  at 
their  separate  positions  (Hooper,  1970).  For  example,  for  an  azimuthal 
density  fluctuation  of  the  form  n(r)  e^m9  where  the  radial  wave  number, 
k^,  is  less  than  the  inverse  probe  separation  (i.e.  kr<[ r  -r  | _1 )  there  is  the 
possibility  of  a  radial  dependence  to  the  azimuthal  phase  shift  as  shown 
pictorially  in  Figure  7  (see  for  example,  Hooper  et  al,  1983  which  is  similar 
to  that  illustrated  here).  With  our  sensor  orientation  along  the  radial 
direction  as  depicted  and  with  the  radial  oscillation  amplitude  of  the  same 
order  of  magnitude  as  the  azimuthal  wavelength  (Ellis  and  Motley,  1974;  Chen, 
1967)  there  is  pictured  a  monotonically  increasing  phase  shift  as  a  function 
of  radius  at  the  two  positions.  Of  course  without  the  addition  of  a  3rd  probe 
in  the  azimuthal  direction  it  is  impossible  to  uniquely  associate  a  propaga¬ 
tion  direction,  k,  with  a  measured  phase  delay.  However,  based  on  physical 
arguments  below,  we  feel  that  Identification  of  our  measurements  with  some 
form  of  drift  mode  propagation  is  reasonable. 

IV.  DISCUSSION  AND  INTERPRETATION  OF  RESULTS 
It  is  well  known  that  a  magnetized  plasma  with  density  and/or  temperature 
gradients  can  support  plasma  oscillations  perpendicular  to  B  with  phase  velo¬ 
cities  near  the  diamagnetic  drift  velocity  (Krall  and  Trivelpiece,  1973). 

For  a  low-8  plasma  with  gradients  in  electron  density  and  temperature 
along  -x  in  a  magnetic  field  B  =  B  z  (see  Figure  6),  the  pressure  balance 


requirement  in  the  absence  of  an  equilibrium  electric  field  yields  a 
zero-order  electron  diamagnetic  drift  of. 
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Plotted  in  Figure  2k  and  2B  are  equilibrium  electron  density  and  temperature 
profiles  along  with  exponential  fits  to  the  data.  From  these  results  and  the 
slab  model  assumption  as  shown  in  Figure  6 
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so  that  the  direction  of  is  along  y  |  7n  |  and  |  7Tg  | 

“n  T 
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are  the  inverse  density  and  temperature  gradient  scale  lengths,  respectively, 
and  we  use  gaussian  cgs  units  with  e  defined  as  the  absolute  value  of  electron 
charge,  K  the  Boltzmann  constant  and  Tg  the  electron  temperature.  If  in  the 
geometry  we  consider,  electrons  are  allowed  to  stream  along  B  (kz  f  0)  in  an 
electrostatic  approximation,  the  Boltzmann  relation  applies  to  each  field 
line  and  is, 
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where  5n  is  the  perturbation  in  denstiy  and  5$  is  the  induced  perturbation 
potential.  Associated  with  the  electrostatic  oscillations  is  an  azimuthal 
fluctuating  electric  field 
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where  we  have  used  equation  (6)  with  the  assumption  =  tQe 
and  Trivelpiece,  1973*  Kelly,  1982).  In  cylindrical  geometry  with  Vng  along 
r,  pressure  forces  are  radial,  drift  wave  propagation  to  lowest  order  is 
azimuthal  and  surfaces  of  constant  density  have  a  twisted  "fluted"  appearance 
(Chen,  1977) .  In  this  description  we  have  implicitly  maintained  the  slab 
geometry  x,y,z  introduced  in  equation  (4)  and  shown  in  Figure  6,  thus  ignoring 
the  cylindrical  curvature  of  the  chamber  to  first  order.  Figure  7  is  a  schematic 
enlargement  of  a  section  of  the  flute  mode  azimuthal  propagation  (Chen,  1977) 
appropriate  to  the  model  at  hand.  Note  that  the  azimuthal  wave  vector 

A 

E  a  kyY  in  the  direction  of  drift  wave  propagation  and  there  is  no  radial 
component  in  this  approximation. 

To  this  point  we  have  sketched  a  zero-order  model  of  drift  waves.  A 
brief  physical  description  of  this  mode  in  relation  to  our  measurement 
technique  is  perhaps  appropriate.  The  nature  of  the  propagation  can  be 
described  as  both  longitudinal  and  transverse  (Chen,  1964);  that  is,  since  it 
and  «E  above  are  collinear  it  is  longitudinal  and,  since  fluid  velocity,  6V 
produced  by  is  perpendicular  to  k  it  is  also  transverse.  Using  the  ideas 
inherent  in  this  model,  6n  and  6<*  are  exactly  y  out  of  phase;  fluctuations  are 
strictly  oscillatory.  However,  under  more  realistic  conditions  (i.e., 
considering  ion  inertia,  polarization  drift  etc.)  the  potential  oscillation  is 
observed  to  lag  the  density  oscillation  in  many  cases  by  an  angle  near  ■jp 
(Hendel  and  Chu,  1970).  In  this  case  wave  growth  will  occur  because,  in  one 
description,  «7  is  outward  at  the  same  time  that  the  plasma  has  already  moved 
outward  (Chen,  1974).  In  the  linear  collisionless  case  this  can  occur  because 
of  the  interaction  between  resonant  electrons  and  the  wave  while  for  the 
collisional  case  it  occurs  because  of  electron-ion  collisions  (Ellis  and 
Motley,  1974).  In  the  linear  treatment,  the  drift  wave  tends  to  be  stabilized 


by  ion  Landau  damping  and  in  the  collisional  regime  by  ion-ion  collisions  and 
the  effect  of  a  finite  Larmor  radius.  Assuming  a  saturation  or  stabilization 
mechanism  for  the  present,  the  relative  potential  fluctuations  can  be 
estimated  from  the  measured  density  variations  as  assumed  in  Equation  (6), 

(i.e.  =  5n  *  It?  (Hendel  et  al»  1968). 

Notice  then  that  for  a  relative  density  fluctuation  of  30*  and  a  thermal 

temperature  of  5000°  (.43  ev),  s  130  millivolts. 

For  the  experimental  arrangement  we  can  consider  the  implications  of 

azimuthal  mode  propagation  in  relation  to  our  measurements.  Along  with 
measurements  of  electron  density  and  temperature  in  Figures  2A  and  2B  are 
exponential  fits  to  each  data  set.  With  the  exponential  dependence  to  these 
profiles  the  temperature  and  density  gradient  scale  lengths  are  constant. 
Therefore  the  electron  diamagnetic  drift  in  Equation  (4)  has  a  radial 
dependence  which  appears  through  temperature  only  (Bv  constant  throughout  the 
measurement  domain).  In  Figure  2C  is  a  plot  of  the  radial  dependence  of 
plasma  potential  as  determined  from  conventional  analysis  procedures  of  the 
E-probe  I-V  characteristic.  Because  of  the  radial  dependence  of  plasma 
potential  there  is  associated  a  radially  outward  electrostatic  electric  field, 
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We  therefore  must  include  any  significant  Doppler  shift  produced  in  the 
laboratory  frame  which  is  associated  with  the  E  x  B  plasma  rotation.  The 


rotation  direction  inferred  from  the  radial  field  is  counter-clockwise  in  the 
cylindrical  frame  and  in  the  ion  diamagnetic  drift  direction  (Hendel  and  Chu, 


1968,  Ellis  and  Motley,  1 97 4 ) .  Hence  a  measurement  of  phase  velocity  in  the 


laboratory,  V  obg,  must  be  corrected  by  the  Doppler  velocity,  7^,  to  determine 
wave  velocity  in  the  plasma  frame  or. 


V  y 

Slab 


(9) 


With  the  assumption  of  azimuthal  mode  propagation  the  observed  velocity,  voba» 
is  related  to  the  observed  spectral  frequency  by, 

I  V  |  =  Sbs  =  “obs  r,  k  =  m 

ky  ra  y  r  (10) 

where  wQbg  is  the  observed  radian  frequency  and  we  consider  an  m  =  1 
rotational  mode.  In  Figure  8  is  plotted  the  frequency  component  of  maximum 
power  for  the  density  fluctuations  as  determined  from  observed  spectral 
measurements.  Although  there  is  some  variation  in  the  observed  laboratory 
frequency,  an  approximately  constant  value  of  40  Hz  is  representative  of  this 
data  set  within  the  measurement  accuracy.  Plotted  in  Figure  9  is  drift  wave 
phase  velocity,  |V  ^1,  in  the  plasma  frame  based  on  the  Doppler  shift  as  seen 
in  Equation  (9)  where  7g.  is  determined  from  plasma  potential  measurements  as 
discussed  above.  Also  shown  is  |7q|,  the  electron  diamagnetic  drift  velocity, 
which  is  determined  through  Equation  (4)  and  the  observed  temperature  and 
density  profiles  of  Figures  2A  and  2B.  From  Figure  9  and  our  assumed  m  =  1 
rotational  model  we  see  that  V^h  *  V^/2  a  result  predicted  by  linear  theory 
for  maximum  growth  of  the  drift  instability  (Hendel  and  Chu,  1968).  Based 
therefore  on  frequency  and  plasma  potential  measurements  and  our  model  of 
azimuthal  wave  propagation  these  results  are  consistent  with  a  drift 
instability  interpretation. 


Recalling  the  model  of  non-local  wave  propagation  discussed  in  relation 
to  Figure  7,  we  now  consider  the  time  delay  associated  with  phase  front  motion 
across  the  position  of  our  two  sensor  system.  This  delay  is  shown  in  Figure 
10  as  a  function  of  radius.  It  is  clear  that  under  an  assumption  of  coherent 
non-local  wave  propagation,  the  measured  time  deiay  is  simply  the  difference 
in  times  of  wavefront  arrival  at  positions  E  and  I.  Under  an  azimuthal  mode 
assumption,  such  as  above,  a  surface  of  constant  phase  necessarily  rotates  at 
a  constant  frequency  in  order  that  the  mode  structure  remain  intact. 

(Notice  that  this  requirement  is  fulfilled  through  the  mode  assumption 
expressed  in  equation  (10)).  Stating  this  somewhat  differently,  the  period 
of  azimuthal  rotation  is  constant  and  independent  of  radius. 
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where  S  (r)  ‘is  the  circumference  and  V  .  s  (r)  the  observed  rotational  velocity 
at  radius  r.  A  measured  time  delay  AT  is  then  the  time  required  to  traverse 
the  distance  AS  in  Figure  7.  (Note  that  this  delay  is  due  to  the  curvature  of 
the  phase  front;  for  rigidly  rotating  "spokes"  the  time  delay  is  zero). 

Because  of  the  constant  rotational  period  the  arc  length  AS  may  be  expressed 
c  (at  (r))  S  (r)  (12) 


Plotted  in  Figure  10  are  arc  lengths,  AS,  based  on  the  velocity 

inferred  from  frequency  measurement  and  the  time  delays,  aT,  observed  at  the 

two  probe  locations  as  a  function  of  radius,  i.e., 

AT  (r)  =  AS  (r) 


Other  than  the  fact  that  these  arc  lengths  are  reasonable  under  model 


(13) 


assumptions,  it  is  difficult  to  compare  the  measurements  to  theory  in  a  more 
complete  manner.  This  is  simply  because  the  model  is  crucially  dependent  upon 
the  curvature  of  the  non-local  wavefront  which  is  assumed  responsible  for  the 
measured  time  delays  and  this  curvature  is  unknown.  It  is  worth  calling  the 
reader’s  attention  further  to  the  fact  that  the  time  delay  and  therefore  the 
curvature  is  not  constant  with  radius;  somewhat  as  if  there  is  an  increasing 
viscosity  or  shear  (Perkins  and  Jassby,  1971)  in  this  view. 

Combining  our  plausibility  arguments  for  a  drift  wave  interpretation  we 
now  consider  the  observed  spectral  distribution  of  ELF  wave  energy  that 
corresponds  to  the  density  fluctuations  which  have  been  the  focus  of  our 
investigation.  Typical  spectral  behavior,  at  different  positions  relative  to 
the  beam  core,  is  shown  in  the  lower  panels  of  Figure  2.  The  observed 
spectral  indices  n  in  an  assumed  power  law  dependence  P  «  f”n  were  typically 
in  the  range  4.5  i  0.3... a  spectral  behavior  consistent  with  drift  wave 
characteristics  (Chen,  1965;  Gary,  1983). 

Finally  we  note  that  the  strongly  magnetized  plasma  electrons  which 

comprise  a  two  component  distribution  (Szuszczewicz,  1983)  of  suprathermals 

and  thermals  in  the  beam  core  and  near  the  beam  boundary,  were  composed 

primarily  of  thermal  particles  in  our  spatial  measurement  domain  outside  the 

beam  core.  In  addition,  we  point  out  that  although  our  zero-order  model  above 

describes  conditions  appropriate  to  the  universal  drift  instability  (kz 

(*0)<<k  ),  an  ion  Larmor  radius  r.  .  of  1.3  u  (T  *  1500°K),  sharp  density 
y  Li  i 

gradients  (L  s  1.4  m)  and  the  introduction  of  collisions  (ion-ion)  suggest 
that  an  assumption  of  magnetized  ions,  as  required  for  positive  growth  of  the 
collisionless  universal  instability,  is  not  rigorously  maintained.  On  the 
other  hand,  the  lower  hybrid  drift  and  ion  cyclotron  instabilities  (Huba, 
Ossakow,  1981a;  Gary  and  Sanderson,  1978)  require,  abiabatic  ions  and 


perpendicular  propagation  (ic  .  B  =  0,  kz  =  0)  for  maximum  growth  so  that  the 

azimuthal  diamagetic  drift  current  is  carried  by  the  ions.  These 

considerations  will  affect  the  conclusions  of  the  equilibrium  model  above  in 

that  the  sign  of  the  drift  wave  frequency  will  change  along  with  Doppler  shift 

considerations  covered  earlier.  Since  the  ion  and  electron  diamagnetic  drifts 

will  be  in  opposite  directions,  and  drift  wave  propagation  is  in  the  drift 

direction,  it  is  straight-forward  to  separate  these  two  possibilities 

experimentally  with  the  addition  of  a  third  probe  in  the  azimuthal  direction. 

In  the  present  series  this  configuration  was  not  in  place  so  that  we  are 

unable  to  unambiguously  determine  a  drift  wave  propagation  direction. 

Under  this  constraint,  however,  and  based  upon  further  characterization 

of  the  plasma,  there  are  some  conclusions  with  regard  to  possible 

instabilities  present.  For  the  particular  series  considered  the  neutral  (N^) 

background  density  can  be  taken  as  6  x  lO1^  cm“^  (P  s  10-^  Torn).  In  our 

temperature  measurement  regime  outside  the  beam  core  (.lev  £  Tg  £  .5ev) 

electron-ion  collisions  dominate  over  electron-neutral  collisions  with 

m  .  s  (5-30)  v  (Braginski,  1965;  Banks,  1966a).  Ion  neutral  collisions,  on 
el  en 

the  other  hand,  dominate  ion  collisions  with  s  (15-130)  These 

collision  estimates  then  provide  some  basis  for  comparison  to  the  collisional 

lower  hybrid  drift  instability.  With  strongly  magnetized  electrons 

(r.  1  cm  at  T  s  .5ev)  and  the  possibility  of  ion  resonance  with  the  wave 

(e.g.  ion-ion  collisions  can  effectively  demagnetize  ions  so  they  are  able  to 

participate  in  resonance  with  a  perpendicularly  propagating  wave),  wave  growth 

can  occur.  For  this  case  one  finds  that  electron  collisions  (v  *  (.2-1)kHz) 

e 

provide  viscosity  to  absorb  wave  energy  and  effectively  place  a  threshold 

condition  on  wave  growth  (Huba,  Ossakow;  1981b).  The  condition  is, 

1/2 


kV.  >  v 
1  e 


where  is  the  ion  thermal  velocity,  ve  is  the  total  electron  collision 
frequency  and  L  is  the  density  gradient  scale  length.  Since  the  model  we 
employ  is  based  on  coherent  non-local  wave  propagation  (L  <  k  ')  and  we  assume 
that  azimuthal  wavelengths  are  comparable  (Chen,  1967),  we  have, 
kVt  ✓  70  sec*1  (k  ^  .001  cm"1,  T1  =  1500°K) 

Therefore  within  our  temperature  range  this  threshold  is  not  reached.  In 
addition,  the  fact  that  this  instability  is  primarily  applicable  in  the  short 
wavelength  regime  (r^^>>.k*1)  would  tend  to  strengthen  an  argument  against  it 
in  the  present  case  where  kr^  ^  0.1.  More  detailed  arguments  must  await  a 
numerical  solution  of  the  dispersion  relation  under  the  conditions  existing  in 
the  chamber  during  the  BPD  measurement. 

Finally,  we  point  out  that  the  instability  destabilization  (or  wave 
growth)  which  occurs  in  transition  to  BPD  can  derive  its  energy  from  a  number 
of  free  energy  reservoirs  available  (i.e.,  temperature  and  density  in 
homogeneity  (presure  gradient),  radial  velocity  shear).  For  our  experimental 
conditions  (krLe  «1,  krLi  $  1,  Ve  >>  the  pressure  gradient  (7ne,  vTg)  is 
perhaps  the  most  likely  candidate  to  contribute  to  destabilization.  For 
example,  in  the  case  of  the  temperature  gradient  the  growth  rate  for  long 
wavelength  collisional  drift  waves  under  the  conditions  above  has  been  shown 
to  vary  with  temperature  gradient  in  Q-machines  (Yamada  and  Hendel,  1978). 


V.  CONCLUSIONS 


We  have  obtained  spectra  of  density  oscillations  in  a  cylindrical 
magnetized  plasma  with  temperature  and  density  gradients  during  the  BPD 
instability  which  occurred  at  an  electron  beam  voltage  of  1.6  kev  and  a  beam 
current  of  40  ma.  We  have  measured  the  electron  density,  temperature,  plasma 
potential  and  cross-correlation  time  delays  associated  with  the  oscillations 
at  separate  radial  positions  using  the  P^  system.  A  zero  order  equilibrium 
drift  wave  model  along  with  an  £  x  B  drift  has  been  compared  to  the  results  of 
measurement.  Based  on  an  m  =  1  mode  assumption,  drift  wave  phase  velocity  is 
estimated  at  near  one  half  the  electron  diamagnetic  drift  speed  if  Doppler 
corrected  velocities  are  used.  The  model  comparison  is  offered  as  a  plausible 
explanation  of  the  observations  based  on  comparisons  to  theory. 
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1  (A)  Side  view  of  the  experimental  configuration  in  the  large  chamber  at 

Johnson  Space  Center  and  (B)  top  view  of  probe  traversal  axis  showing  loca- 
i  with  respect  to  the  electron  beam  gun  and  target. 


—  Electron  (Igg)  and  Ion  (Iqj )  saturation  currents  collected  by  the  E . 
respectively.  Ion  currents  have  been  multiplied  by  -1  for  comparison 
ise  difference  between  Idr,  Iri. 
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Fig.  4  —  Cross-correlation  integral,  C(r),  of  normalized  electron  and  ion  saturation 
current  variations  using  the  IgE>  *BI  segment  shown  in  Figure  2. 
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Fig.  8  —  Major  frequency  component  of  spectra  as  a  function  of 
radial  position  with  m=l  calculated  results 
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Fig.  9  —  Measured  electron  diamagnetic  drift  VD  and  measured  Vph  (corrected  for  Doppler) 

radial  position 
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Fig.  10  —  Measured  time  delay  between  the  two  probe  positons  from  cross  correlation  results 
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